Orbital Switching and the First-Order Insulator-Metal Transition in 

Paramagnetic V2O3 
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The first-order metal-insulator transition (MIT) in paramagnetic V2O3 is studied within the ab-initio 
scheme LDA+DMFT, which merges the local density approximation (LDA) with dynamical mean 
field theory (DMFT). With a fixed value of the Coulomb U = 6.0 eV, we show how the abrupt 
pressure driven MIT is understood in a new picture: pressure-induced decrease of the trigonal 
distortion within the strong correlation scenario (which is not obtained within LDA) . We find good 
quantitative agreement with (i) switch of the orbital occupation of (oi g , eL, e^ 2 ) and the spin state 
S = 1 across the MIT, (ii) thermodynamics and dc resistivity, and (Hi) the one-electron spectral 
function, within this new scenario. 
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PACS numbers: 71.28+d,71.30+h,72.10-d 

The spectacular metal-insulator transition (MIT) in 
the paramagnetic (P) state in V2O3 has long been con- 
sidered as a classic, and by now almost a textbook ver- 
sion of correlation-driven MIT in a one-band correlated 
system described by the one-band Hubbard model |Q. 
This conventional wisdom has been recently challenged 
by new experiments, clearly showing the coupled spin- 
orbital character of the system, and necessitating a revi- 
sion in terms of a multiband picture. 

The conventionally accepted picture rested upon the 
following argument. In the high-T phase, V2O3 exists in 
the corundum structure, and with a 3d 2 state (V 3+ ), the 
two orbitals are empty, while the triply degenerate t2 g 
orbitals are filled by two electrons. This triple degeneracy 
is lifted by a small trigonal distortion, leading to a singly 
occupied a\ g orbital oriented along the c-axis, and dou- 
bly degenerate planar orbitals, occupied by the sec- 
ond electron. Castellani et al. || proposed that strong 
covalent effects lead to a bonding singlet a\ g state in- 
volving two V ions along the c-axis. The remaining elec- 
tron in the two-fold degenerate orbitals gives rise to 
a S = 1/2 model with orbital degeneracy. This inspired 
development of theoretical techniques, culminating in the 
dynamical mean field theory (DMFT) leading to con- 
siderable improvement in our understanding of the MIT. 
Within the basic picture (2), Rozenberg et al. ||] used 
DMFT for one- and two-orbital Hubbard models with 
Bethe density-of-states to study the MIT in V2O3. 

Recent polarized X-ray scattering results of Park et 
al. [^| require, however, an interpretation in terms of 
a spin S = 1 at each V site, with a mixed orbital 
e^a\g : e K g e g ~ — x : (1 — x) configuration. An exciting con- 
clusion from these results is that the above ratio changes 
its value abruptly at the MIT, forcing one to abandon 
the one-band Hubbard model to describe V2O3. Notice 
that this implies an important role for the trigonal split- 
ting, since the lower-lying orbital will be more "localized" 
when local Coulomb interactions are switched on. This 
raises questions concerning a possible link between the or- 



bital "switching" and the drastic change in the electronic 
state, and to a possible common underlying origin. 

The high-spin ground state results from strong Hund 
coupling, a fact borne out by LDA+U ||. Furthermore, 
LDA+U results argue in favor of models without orbital 
degeneracy. This line of thought has been extended in 
detail by Tanaka jf) , including spin-orbit coupling effects 
in a cluster-based approach. Mila et al. JsJ have pro- 
posed a model with 5 = 1 and orbital degeneracy for the 
antiferromagnetic insulating phase. 

Hallmarks of strong electronic correlations in V2O3 
have been observed across the MIT in thermodynamic [|| , 
dc transport |l(J , photoemission |ll| , and optical [jl2| re- 
sponses. In particular, specific heat measurements give 
7 ei ~ 4.4, p dc (T) = pa + AT 2 with large A in the PM 
phase, while PES and optical measurements clearly re- 
veal gap formation in the PI phase, and the development 
of low-energy quasicoherent weight in the PM, with a 
dynamical spectral weight transfer (SWT) from high- to 
low energy over a wide energy scale (~ 4 eV). This type 
of dynamical SWT-driven MIT is beyond the scope of 
LDA+U, which treats dynamical effects of local Coulomb 
interactions within static Hartree-Fock (HF) approxima- 
tion. These are precisely the effects captured reliably by 
DMFT. 

Thus, a consistent description of the MIT, along with 
an understanding of the strong correlation features re- 
quires a combination of structural aspects of V2O3 (en- 
coded in LDA+U) with a reliable many-body theory like 
DMFT. In this letter, we study these questions within 
LDA+DMFT which has been shown to provide a 
good quantitative ab-initio description of correlated elec- 
tronic systems. Moreover the MIT in the P-phase has 
been studied |J] using LDA+DMFT(QMC). However, 
the link between the MIT and the abrupt switching of 
orbital occupation has, to our knowledge, not been ex- 
plored in detail. We should mention that such a scenario 
may have broader application to other systems, notably 
in Ca2- x Sr x RuO± fl4[], and is a hallmark of the 
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FIG. 1. LDA partial density of states for the (dashed) 
and dig (solid) orbitals, obtained from Ref. [13]. 

importance of orbital correlations in a system. 

We start with the actual LDA bandstructure of V2O3 
in the corundum structure jDJ. These results repre- 
sent (Fig. metallic behavior, with comparable width 
(2.5 eV) for both ai g and e g bands, and a strongly asym- 
metric structure. With an onsite U ~ 5 — 6 eV, this 
would invalidate a MO-based approach || . Further from 
the aig — splitting, the trigonal field splitting is es- 
timated to be ~ 0.4 eV. The LDA densities of states 
with and without the trigonal distortion do not differ 
much, and in fact are both representative of good metal- 
lic behavior, showing that the MIT cannot be related 
purely to A trg . Notice that at the local level, external 
pressure decreases A trg , resulting in a sudden increase 
in n aig around A trg — 0.28 eV 0. Turning on the lo- 
cal interactions, U, Jh and U' = (U — 2Jh) will push up 
most of the spectral weight to the Hubbard sub-bands, 
and strong Jh will favor S = 1 at each V site. So, 
while these "ground state" features are reasonably well 
accounted for by LDA(+U), aspects like the mixed or- 
bital configuration, the switch in orbital occupation at 
the MIT accompanied by dynamical SWT require a reli- 
able treatment of dynamics of correlated electrons. 

Thus, the one-particle part (LDA) of the Hamiltonian 
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with the total DOS, p° total {e) = Eka/3 CT 8(e-e a/3 (k)), and 
a,f3 = a\ g , e^, e g2 . To avoid double counting of interac- 
tions which are already treated on the average by LDA, 
we follow to write 
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where 



U(n aa - 5) + \Jii{ri aS - 1), with 



U, Jh as defined below. 

With the interactions in the ti g sector, the full Hamil- 
tonian reads, 



H = Hn 



u 



(3) 



Constrained LDA calculations yield U — 5 — 6 eV (with- 
out inclusion of screening of t2 g interactions by e g elec- 
trons), J H ^ 1 eV, and U£%' = U' = (U - 2J H ) = 
3 — 4 eV. Following we use the fact that the e a g 
bands are well separated from the ti g bands to consider 
only the t2 g manifold. Furthermore, in the PM and para- 
orbital phase, we have G a p aa > (u>) = <5a/3<W'G Qcr (cL>) and 

In the t2 g sub-basis, a DMFT solution involves (?) re- 
placing the lattice model by a self-consistcntly embedded 
multi-orbital, asymmetric Anderson impurity model, and 
(ii) a selfconsistency condition which requires the im- 
purity propagator to be equal to the local (/c-averaged) 
Green function of the lattice, given by 
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Using the locality of T, a p in d — 00, we have G a (uj) = 
G° (u> — E a (w)) from the Hilbert transform of the LDA 
DOS. Also, importantly, the inter-orbital couplings scat- 
ter electrons between the a\ g , e^ bands, so that only the 
total number, nt 2g = J2 a nt 2 g ,a i s conserved in a way 
consistent with Luttinger's theorem. 

To solve the multi-orbital, asymmetric Anderson im- 
purity problem, we use the iterated perturbation theory 
(IPT), suitably generalized to the case of t2 g orbitals for 
arbitrary filling fll5| . The local propagators are given by 
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Local self-energies S Q (w) are computed within an ex- 
tended IPT scheme that explicitly satisfies the gener- 
alized Friedel sum rule (Luttinger's theorem) to a very 
good accuracy. Mathematically, 
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where, for example, 
U 2 

= N ai -^- G° a (iw„)G°(iw m )G°(»w B + iuj m - iw) 
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with N a ~ ( — 2 for a, 7 = e^e^ and 4 for a, 7 
ai g ,Bg 12 . The bath propagator is G a {oj) = [oj + (x a 
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numbers defined from G a and G° a , and the inter-orbital 
correlation function is D al [n] — {n a n 7 } — (n a )(n 7 ) — 
f+™ f(w)Im\E a (u)G a {u)]du. This last identity 
follows directly from the equation of motion for G a (uj). 

These coupled equations are solved selfconsistently to 
obtain the spectral function. We choose A trg — 0.32 eV, 
completely consistent with the LDA. To study the MIT, 
we notice that external pressure decreases the quantity 
Atrg = E aig — E e 7r and leads to an increase in n aig ; in 
particular, it leads to a reduction in A trg across 0.28, at 
which point, a sudden increase in n aig has been reported 
from cluster calculations 0. To study this effect, we 
monitor the spectral function for different values of n ai , 
with fixed total number of electrons. 

In Fig. ||, we show the partial DOS for the e" g ,ai g 
orbitals, for the metallic (dashed) and "Mott" insulat- 
ing (solid) phases, corresponding to total a± g occupation 
n aig = 0.41, 0.36, respectively. As expected from the trig- 
onal splitting assignment, the lower- lying orbitals are 
more localized in the solid, with a gap, A e * = 0.45 eV. 
Correspondingly, the a\ g band is more itinerant, com- 
pletely consistent with the fact that t ai is by far the 
largest hopping integral in the real system. From Fig. ||, 
we estimate Ai ns = 0.35 eV. With n a ,„ = 0.41, the 



while the a\ g spectrum develops a very narrow, quasi- 
coherent peak with a FWHM = 0.07 eV at E F . We 
interpret these results as a microscopic derivation of the 
phenomenological "two-fluid" models proposed earlier to 
understand metal- insulator transitions ||. In our new 
picture, increasing pressure decreases A trg , resulting in 
increased population of the a\ g orbital, and leading (via 
ta lg ,a lg ) to an abrupt MIT accompanied by large dyna- 
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FIG. 2. t 2g partial density of states using LDA+DMFT for 
U = 6.0 eV in the PI (solid) and PM (dashed) phases of V 2 3 . 
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FIG. 3. The integrated inverse- and direct photoemission 
lineshapes for the PI (solid) and PM (dashed) phases within 
LDA+DMFT for U = 6.0 eV. The LDA result is the thin 
solid line. 

mical SWT. 

Using the corresponding self-energy, we carry out a 
quantitative estimate for the effective mass enhancement 
and the yl-coefficient of the low-T quadratic term in the 
dc resistivity. Indeed, from i?e[E 0l9 (w)], we estimate 



spectral function still exhibits insulator-like features, m*/m — [1 



4.16, close to the value of 4.4 



extracted from low-T specific heat measurements. Also, 
A = Pdc (T)/T 2 = (m*/ne 2 )\d 2 ^ lg (u)/dLu 2 \ is large, be- 
cause \d 2 Y!' aig (uj) I duj 2 \ = 38 > 1. This will result in a 
large quadratic term in the low-T dc resistivity, as ob- 
served. 

In Fig. we show the integrated PES and IPES spec- 
tra in the metallic and insulating phases, along with the 
corresponding LDA spectra. As one would expect, and 
consistent with observations, the PES spectrum shows a 
very narrow quasicoherent peak, with most of the spec- 
tral weight in the incoherent "lower Hubbard band" , cen- 
tered at u> = —3.0 eV. The PES intensity at Ep is 
strongly reduced from the LDA prediction. Quantitative 
comparison with the PES spectra measured at higher T 
needs an extension of our approach to finite T; however, 
we expect that the strong T-dependence of the quasi- 
coherent feature within DMFT will smear out the low- 
energy peak for T > T co h — 300K, bringing the spectra 
in closer agreement with published PES results. Signa- 
tures of strong correlations are also visible in the IPES 
spectra as "upper Hubbard band" features, and the dif- 
ference between the LDA and LDA+DMFT spectra is 
striking. 

Especially interesting is the change in the spectral 
function at Ep as a function of n aig . In Fig. |], we show 
this variation for our chosen parameter set. Astonish- 
ingly, p(Ep) exhibits a sharp jump from to 0.261 
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FIG. 4. The total DOS at Ef as a function of the oi s 
orbital occupation. 

around n Qls = 0.39, clearly showing that the MIT is in- 
deed first order, in agreement with observations. Our 
results imply that the MIT is accompanied by an abrupt 
jump in the itinerant carrier density, rather than increase 
in carrier mobility, and is consistent with the conductiv- 
ity jump at the PI-PM transition. From the computed 
spectral functions, (n aig ,n e * i ,n e * 2 ) — (0.36,0.82,0.82) 
in the PI, and (0.41, 0.79, 6.79) in the PM. Our cal- 
culated orbital occupations are slightly different from 
those estimated by polarized XAS measurements || , but 
show the correct trend across the MIT. Further, we find 
( s fz) = zJ2 a ([ n ial ~ n ia [\ 2 ) = 0.92 for the chosen pa- 
rameter set, consistent with the observation of S = 1 at 
each V site. Finally, we monitor D aig ^ and D e ^^^^ 
and find that both decrease across the MIT: in fact, 



D aig 



= -0.043 for 



= 0.36 (PI), and equals 



-0.007 for n aig = 0.41 (PM), while ^ = -0.065 
(PI) and -0.036 (PM) showing clearly that the MIT is 
accompanied by a strong reduction of local (anti-ferro) 
orbital correlations. 

As mentioned earlier, the MIT in V2O3 has recently 
been studied ||j by LDA+DMFT(QMC), so it is in- 
structive to compare our approach with theirs. Though 
QMC is numerically exact, it cannot access the low- 
T regime (T < 400K). Further, analytic continuation 
of Matsubara (imaginary) frequencies onto the real fre- 
quencies is a highly ill-posed numerical problem. While 
the IPT is not exact, it is free from these drawbacks, 
and is known to yield reliable results for the problem 
at hand fl5|| , enabling one to study thermodynamic and 
transport properties at relatively modest computational 
cost. In contrast to the earlier pioneering study, we 
have explicitly shown the existence of the first-order PI- 
PM transition and its intimate link to the accompany- 
ing orbital switching in V2O3. We have not attempted 



to make a detailed comparison with finite-T > 170 K 
spectral functions in the PI and PM regions, a program 
which has been carried out in Ref. |D| (but only for 
T > 400 K, where the first-order MIT is replaced by a 
smooth crossover), and we plan to report these aspects, 
along with finite-T transport properties and lattice ef- 
fects, in a future work. 

In conclusion, we have presented a different scenario, 
based on combination of LDA bandstructure with dy- 
namical, local spin and orbital correlations, to under- 
stand salient features of the first-order MIT in the P- 
phase of V2O3. The first-order MIT is shown to be accom- 
panied by changes in orbital occupation in a way quali- 
tatively consistent with observations. Further, very good 
agreement with low-T thermodynamic and dc resistivity 
is obtained within the same approach. This represents a 
new picture for the correlation-driven MIT in V2O3; one 
which has applications to other systems where coupled 
spin-orbital correlations result in changes in orbital occu- 
pations across the MIT. Such an approach should be ap- 
plicable, with extensions to include broken symmetries in 
the spin/orbital sectors, to systems like Ca2- x Sr x RuO^ 
as well. 

We are indebted to L. H. Tjeng for many enlightening 
discussions. Work carried out with the support of the 
Sonderforschungsbereich 608 of the Deutsche Forschungs- 
gemeinschaft. 
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